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' Two Body Dirac Equations (TBDE) of Dirac's relativistic constraint dynamics have been suc- 

cessfully applied to obtain a covariant nonperturbative description of QED and QCD bound states, 
r-y . Coulomb-type potentials in these applications lead naively in other approaches to singular relativis- 

tic corrections at short distances that require the introduction of either perturbative treatments or 
C ' smoothing parameters. We examine the corresponding singular structures in the effective potentials 

of the relativistic Schrodinger equation obtained from the Pauli reduction of the TBDE. We find 
that the relativistic Schrodinger equation lead in fact to well-behaved wave function solutions when 
the full potential and couplings of the system are taken into account. The most unusual case is the 
coupled triplet system with S = 1 and L — {(J — 1), (J + 1)}. Without the inclusion of the tensor 
coupling, the effective S'-state potential would become attractively singular. We show how including 
O |' the tensor coupling is essential in order that the wave functions be well-behaved at short distances. 

For example, the S'-state wave function becomes simply proportional to the D-state wave function 
i-C . and dips sharply to zero at the origin, unlike the usual 5-state wave functions. Furthermore, this 

behavior is similar in both QED and QCD, independent of the asymptotic freedom behavior of 
\ the assumed QCD vector potential. Light- and heavy-quark meson states can be described well 

by using a simplified linear-plus-Coulomb-type QCD potential apportioned appropriately between 
■ world scalar and vector potentials. We use this potential to exhibit explicitly the origin of the large 

| n-p splitting and effective chiral symmetry breaking. The TBDE formalism developed here may be 

£ — , used to study quarkonia in quark-gluon plasma environments. 
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PACS numbers: 12.39.Ki,03.65.Pm,12.39.Pn 

I. INTRODUCTION 

The two-body Dirac equations we discuss in this paper are based on Dirac's constraint formalism and a minimal 
interaction structure for two particles in relative motion, first used by Todorov [1] and confirmed by both classical 



[2] and quantum field theory [3]. The constraint approach gives more than a sophisticated method for guessing 
relativistic wave equations for systems of bound quarks or general fermion-anti-fermion systems, since it can be 
readily combined with the field-theoretic machinery of the Bethe-Salpeter equation. When used with the kernel of 
the Bethe-Salpeter equation for QED, it combines weak-potential agreement in QED [4] with the nonperturbative 
structure of the field-theoretic eikonal approximation [1, 3]. The minimal interaction structure is then automatically 
inherited from relativistic classical [5] and quantum field theory [3]. 

As has been demonstrated earlier, the constraint equations correspond to a "quantum-mechanical transform" [6, 7] 
of the Bethe-Salpeter equation (BSE). This is provided by the two coupled Dirac equations whose fully covariant inter- 
actions are determined by QED in the Feynman gauge [3, 4, 8]. Unlike most other truncations of the Bethe-Salpeter 
equation, the constraint approach does not require the use of the awkward Coulomb gauge (whose noncovariant nature 
does not allow its incorporation in covariant equations). Instead, its expansion about the BSE naturally occurs in 
the covariant Feynman gauge and is free from spurious infrared singularities that occur in the other approaches when 
that gauge is used [9]. 

In QCD with flavor-independent interactions this formalism leads to spectral results in very good agreement with 
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most of the experimental meson spectra 1 . At the same time, as stressed in a recent publication [10, 11], the formalism 
naturally accounts for the perturbative results of QED bound states, when treated in a nonperturbative manner. So 
far this has not been fully replicated in any other approach. In a natural way it leads not only to good singlet-triplet 
state splittings for the light as well as heavy mesons, but also to a Goldstone behavior for the pion. By this we mean 
that the numerically computed pion mass tends to zero when the quark mass tends to zero. This is tied to the same 
relativistic structures that account for the nonperturbative positronium and muonium results [10]. 

The relativistic Two-Body Dirac equations may be written as an effective one-body wave equation [12-14]. The 
proper formulation of this relativistic scheme requires the successful treatment (that is a covariant elimination) of 
the quantum ghost states (due to the presence of the "relative time") that first appeared in Nakanishi's work on the 
Bethe-Salpeter equation [15]. These coupled constraint equations, known as the Two Body Dirac Equations (TBDE), 
have the following important characteristics. Firstly, in the special limit in which one of the particle masses becomes 
infinite, the equations reduce to the (one-body) Dirac equation. Secondly, in the general case the Pauli reduction of 
the TBDE leads to a relativistic Schrodinger equation, which is the same as the non-relativistic Schrodinger equation 
in form but is relativistic in content, including all relativistic spin and relativistic kinematics. It displays various 
spin-spin, spin-orbit, tensor and Darwin terms with energy-dependent denominators. Relativistic kinematics have 
been properly taken into account to give a relation between the eigenvalue of the relativistic Schrodinger equation 
and the invariant mass of the composite system. 

We explore several related questions in this paper. How is it that the TBDE or its equivalent relativistic Schrodinger 
equation leads to a Goldstone-like behavior of the pion? How does this approach bypass many of the singularities that 
appear in the effective potentials in other approaches without the necessity for introducing cutoff parameters? A more 
complete understanding of how the TBDE are able to accomplish this not only will aid in a better understanding of 
its success in meson spectroscopy but will assist us in its application to two-body bound states in other environments 
such as those in a strongly-coupled quark gluon plasma. 

In the application of TBDE to QED and QCD bound state problems, the interaction includes Coulomb-type poten- 
tials which lead naively in other approaches to singular relativistic corrections at short distances (delta functions and 
potentials more attractive than — l/4r 2 ) that require the introduction of either perturbative treatments or smoothing 
parameters. In the weak potential limit in which the potential is regarded as small compared with the masses and 
center of momentum (cm.) energy, the effective potentials of the relativistic Schrodinger equation, obtained from the 
Pauli reduction of the TBDE also displays these types of singular potentials. However, as we shall see these potentials 
become nonsingular when one uses the strong potential form in which one does not ignore the potential compared with 
the masses and cm. energy. Nevertheless, the potential of the triplet system, 2S+1 Lj, with S = 1 and L = (J — 1) 
and ( J + 1), remains much more pernicious as the potential retains its singular behavior independent of whether one 
regards the potential as weak or strong. We would like to describe here these unusual singularity structures in the 
potential and show how the TBDE formalism leads to wave function solutions that are nonetheless well-behaved and 
physically acceptable, when the full couplings of the system are taken into account 2 . 

In addition to our investigation of the singularity structures in TBDE, we also wish to test a simplified QCD 
based potential for future applications. Previous studies of the TBDE [8, 17, 18] made use of the QCD based Aldcr- 
Piran potential [19] giving a quite successful description for various qq states. The Adler-Piran potential contains 
functions with many terms and parameters divided into different sections of the spatial region. On the other hand, 
simple potentials such as the Cornell potential [20] appears to be adequate for many applications, although lacking 
in asymptotic freedom. It is therefore desirable to seek a simplified potential similar to the Cornell potential, but one 
containing asymptotic freedom that can be easily modified for future application of the TBDE in other quarkonium 
problems. One such application is on the stability of qq states in the quark-gluon plasma which consists of quarks 
and antiquarks of different flavors, and gluons. The degree to which the constituents of a QGP can combine to form 
composite entities is an important property of the plasma. 

To facilitate the application of the TBDE in terms of the equivalent relativistic Schrodinger equation, we present 
relevant useful details to indicate how various Darwin and spin-dependent potential terms can be constructed in 
Appendix C. Once the various terms of the potential have been constructed, the solution of bound states problems 
in the TBDE is mathematically just as simple as the solution for bound state problems in non-relativistic quantum 
mechanics. 

Accordingly, we begin in Sec II with a discussion of the most often cited QCD potentials (including a simplified 
model for the QCD potentials to be used for the first time in this paper in conjunction with the TBDE for meson 
spectroscopy). It is well known how a naive use of Coulomb- type potential leads to singularity problems when 



1 Isoscalars such as the r], r)' , and w are not included 

2 We note that the finite- r singularity structures that occur in the Breit equation do not appear in the TBDE [16]. 
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one introduces relativistic spin-spin corrections. How we avoid these problems in the TBDE is the next question 
considered. In Sec. Ill, we seek out first how these problems are avoided in the Dirac equation for an external 
Coulomb potential. This allows us to explore in Sec. IV the parallels between the singularity structures that occur 
in the Dirac and TBDE and how they each connect the different structures of the weak-potential (or perturbative) 
and strong-potential (or nonperturbative) forms of the two sets of equations. Those different structures, although 
requiring different treatments, give the same spectral results. By this we mean that a perturbative treatment of the 
weak-potential forms gives the same spectral results as a nonperturbative (analytic or numerical) treatment of the 
strong-potential or nonperturbative forms of the two sets of the equations. The example we use in Sec. V for this 
parallel discussion is the TBDE for electromagnetic interactions, QED. In Sec. VI we discuss the spectral results, 
focusing on the pion, the p, and the singularity structure of the TBDE for a simplified QCD potential model in the 
case of 3 Si- 3 Di or more generally [ 3 (J — l)j] -[ 3 (J + l)j] mixing. In Sec. VII we summarize the results and discuss 
questions that may arise when attempts are made to apply the TBDE to two-body bound states in a QGP. 



II. QCD MODEL POTENTIALS 



Previously, the authors of [17] used a sophisticated form of the static quark potential developed by Adler and Piran 
[19], one that has ties at all length scales to field theoretic data. Very good agreement with experimental quarkonium 
spectrum was obtained. On the other hand, in nonrclativistic treatments the most commonly used static quark 
potential for potential model studies is the Cornell potential [20], 

V(r)=-^+br, (1) 

as in [21, 22]. Although not displaying asymptotic freedom, it does give the dominant Coulomb-like behavior as well 
as the linear quark confinement. Early on a model was proposed by Richardson for a static potential which a) depends 
only a single scale size A, and b) interpolates in a simple way between asymptotic freedom and linear confinement 
[23]. Richardson's model for the static interquark potential in momentum space is 

~ . . 167T 1 

(q) = _ ~2Tq 2 ln(l+ q 2 /A 2 )' [) 

arising from the assumption that 

3q 

(including the color factor —4/3). Asymptotic freedom requires that for q 2 /A 2 >> 1 



^q)=-^F, (3) 



, 2 \ 127T 1 

as(q) ^^MqW W 



On the other hand, the property of linear confinement requires that for Ar >> 1, V(r) a r or cquivalcntly that 
for q 2 /A 2 << 1 one must impose a s (q 2 ) <~ q . The interpolation of Eq. (2) is not tied at all in the intermediate 
region and only roughly tied in the large r region to any field theoretic data. Nevertheless it provides a convenient 
one-parameter form for the static quark potential. In coordinate space it has the form 

Vir)-*gZ-*l£>, (5) 

where /(Ar) is given by a complicated integral transform 3 that displays the asymptotic freedom behavior for r — > 
of 

/(AO - ~, (6) 

while for r — > oo, 

/(Ar) 1. (7) 



In addition to the spin independent nonrclativistic model presented in [23] see also a relativistic extension of it given in [24]. 
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A simpler model, which we will apply in this paper and one which displays the same large and small r behavior is 4 

_ 87rA 2 r 16tt 
W ~ ~Yt 27rln(e 2 + l/(Ar) 2 )' ( ' 

It amounts to replacing Richardson's /(Ar) by 2/ln(e 2 + l/(Ar) 2 ), having the same limits. Although not giving 
as good a fit to the spectra as the more closely tied QCD based potential of [19], the modified form of Eq. (8) 
which we use in this paper (see Eq. (52) below) does provide reasonable results for the spectrum. Furthermore its 
linear-plus-Coulomb-type parametrization is more convenient for extension of the quark model to high temperature 
environments. 

Problems arise in the quark model with the above potentials if their relativistic corrections are naively grafted from 
semirelativistic expressions. For example, the spin-spin interaction 



V 2 W r er 2 

6mi77l2 



(9) 



would lead to a singular delta function potential that can only be treated perturbatively. Some approaches simply 
include cutoff parameters so that the Laplacian is not singular. How do the TBDE treat this problem? Potential 
energy terms such as the above arise from the second order reductions of those equations (the Pauli forms). Let us 
first examine how such problems are treated in a very natural way in the Schrodinger Pauli-form of Dirac's original 
wave equation. 



III. SINGULARITY STRUCTURE OF PAULI-FORM OF THE DIRAC EQUATION 

The q-q interaction in Eq. (8) contains the color-Coulomb term that is proportional to 1/r and a logarithmic function 
of r. It leads naively to singular relativistic corrections that may render the solution singular at short distances. It is 
worth while to investigate Coulomb-type potential in relativistic equations. Let us be more precise in our definitions 
of singular potentials. Case [26] describes how potentials that are more attractive at the origin than — l/4r 2 must be 
adjusted to maintain their self-adjoint status. Let us call such potentials, attractive singular potentials. They include 
attractive delta functions and attractive 1/r 3 potentials that appear in spin-orbit terms. Such terms must either be 
treated only in perturbation theory or in cases where the coupling is strong, require adjustments, e.g. by smoothing 
parameters. Calogero [27] and Frank et al. [28] also discuss another category of potentials called repulsive singular 
potentials. These are repulsive potentials that exceed an inverse quadratic power law behavior. Strictly speaking they 
need not be treated using perturbation theory, although in the case of weak potentials they are most easily treated 
like that. We will discuss these more in the context below. 

Let us show how the Pauli-form of the Dirac equation with a Coulomb-type interaction contains effective potentials 
that are repulsively singular, when viewed in an incomplete or perturbative context. However, when viewed in a 
complete or nonperturbative context the effective potentials are nonetheless nonsingular. We examine for simplicity 
the case of the Dirac equation in a Coulomb potential (A = —a/r) (instead of the more complicated forms with 
asymptotic freedom) for stationary states, 

(a ■ p+/3m + A)$ = E^S. (10) 

Then with 



we have 



Eliminating \ we obtain 



m — E + A rr p 

erp -E-m + Al U 



0. (12) 



4 An earlier coordinate space form that displays asymptotic freedom as well as linear quark confinement proposed in [25] is V = (87r/27)(l — 
Ar) 2 /(r In Ar). 
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leading to 



2 iA' A' 



( p - TFT ZT f ' P+ TFT T\ (T ■ L ^ = ^ E A r m ^ ( 14 ) 

(E + to — A) r(E + m — A) 



We can eliminate the first order derivative f • p term by the substitution 

4> = F(r)ip (15) 

if one takes 

F' A 1 



F 2(E + m-A)' 

Then our equation becomes 



(16) 



< P ° + HE^m-A) + Z iEV m -A)> + r(E + rn-Af ' ^ = [{E ~ ^ " m ^ ^ 
which for a Coulomb potential becomes 

. 2 2Ea a 2 2ira5 3 (r) 3 a 2 a , 2 2n , , N 

(P - ^- ^ + (^ + m + a/r) + 4 r 4 (g + m + a/r) 2 + r^E + ra + a/r) °" ' = " ™ ^ (18) 

For a > 1/2 the inverse quadratic term would lead to an overall attractive singular potential for 5-states. If one 
takes the weak potential limit in which the denominators in the three succeeding terms are replaced by 2m then the 
spin-orbit potential is an attractive singular potential for both coupling states and must be handled by perturbative 
techniques. Under weak potential circumstances the delta function potential would be treated by perturbative 
techniques. Since it is repulsive, one could, in principle treat it in a nonperturbative way. It has been shown, however 
[29] , that a nonperturbative treatment of repulsive delta functions potentials produce no effect on bound state energies. 
The repulsive 1/r 4 term would require special numerical treatments (its perturbative effects on S-states is ill-defined) . 
Let us now compare perturbative and nonperturbative treatments of this Pauli form. 

By using the atomic units r = x/(Ea) the above equation takes the dimensionless coordinate space form of 

2 2 2o? a 4 2 7 m 4 £ 3 (x) 
k x x x 2 + (l + m/E + a 2 /x) [ ' 

3 a 6 a 4 

+ 4 x 4 (1 + m/E + a 2 /x) 2 + x 3 (l + m/E + a 2 fx) " ' 



V 2 A 3 (VA) 



1)V 



The standard perturbative treatment retains terms through order a 4 to arrive at the equation, 

( - - 2j f-^+ ™ 4 * 3 w+|^ • w = [i - Q V (20) 

The standard semirelativistic spectral results through order a 4 can be obtained by treating this as an ordinary 
eigenvalue problem with the last three terms on the left hand side as a perturbation. For the ground state this leads 
to 

„ ma 2 ma 4 «. 

E = m-— — + 0(a e ). (21) 

Note that the a 6 term in Eq. (19) does not contribute perturbatively to this order. 



A. Small r Effective Potential and Wave Function Behaviors 



The <5 3 (r) potential in Eq. (20) must only be treated perturbatively to obtain a nonzero result. We know, however, 
that the Dirac equation in this case can be solved analytically. How does that reconcile here with the appearance of 
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these singular potentials, particularly the S 3 (r) potential? Let us restrict ourselves here to S-states to make our main 
point. That would mean with u = y/iwxtl) that we must include all terms in 



2 d 2 2a 2 a 4 3 a 6 _ f m \ 2 ] 

{ ~ a d^~~~^ + 4x 2 (x(l+ m /E) + a 2 ) 2 > U - [1 ~{E~) K (22) 

in any nonperturbative solution. Note that we have left out here the <5 3 (x) term in Eq. (19) since it, together with 
the Coulomb potential in the denominator would yield a vanishing result for its contribution (its expectation value 
in any well-behaved basis would give zero). That means that we must have the rather unusual circumstance here 
of the term which does not contribute to the weak potential form Eq. (20) (the a 6 term) having a nonperturbative 
effect on the spectrum that reproduces that of the perturbative (5 3 (x) term (in a perturbative expansion). A set of 
straight-forward but tedious manipulations show how this comes about. We first point out that the a 6 term in Eq. 
(19) has a short distance behavior of that of a repulsive r~ 2 behavior that is lower order in a, so unlike its weak 
potential 1/r 4 form it is not in the category of a repulsive singular potential. This allows a standard type of solution. 
One finds that the ground state wave function and eigenvalue are given by 

u = kx p (x(l + m/E) + a 2 y cxp(-Ax), 
1 

7 = -o, 



E = m\J 1 — a 2 , 

x = J3SE± = i (23) 

a VI - or 

We also verify that our exact solution 

E = mV 7 ! -a 2 =m(\-\ 2 a 2 - ^a 4 + 0(a 6 )^j , (24) 

agrees with the perturbative spectral results Eq. (21). Note that the small r behavior of the radial part of the wave 
function is 

V> ~ x^^- 1 / 2 (25) 

which dips toward the origin unlike the flat behavior of the nonrelativistic limit of the Pauli form or mildly singular 
behavior of the Dirac wave function. 

In summary, the Pauli-form (19) of the Dirac equation in the weak potential approximation or the perturbative form 
of (20) includes terms missing in the strong potential or the nonperturbative form (22 ). Also, Eq. (22) includes terms 
missing in the weak potential or the perturbative form (20). However, they both give rise to the same spectral results 
through order a 4 , with one treated in a weak potential approximation and the other treated with no approximations 
(whose spectra expansion yields the same result). As we shall see below, such an unusual feature (with different 
parts of the equation contributing to the perturbative and nonperturbative spectral evaluations) is also displayed in 
the Pauli-form of the TBDE of constraint dynamics. 

The above exercise shows that what appears as singular in a perturbative context turns out in fact to be non- 
singular in a full non-perturbative treatment. For our case of the ground state, the two approaches give the same 
result up to order a 4 , using different parts of the effective interaction. Does this extend to the radially and orbitally 
excited states? Do different parts (and approximations) of the interaction used in reaching the spectral results lead 
to the same results for all high excited states. Although we have not shown this here, it is expected to be true since 
the exact (nonperturbative) Pauli form should faithfully reproduce the exact spectral results of the first order form 
of the Dirac equation. This would imply that one would expect the two approaches to give the same results through 
order a 4 . 



5 At short distance we have -j^ + ( 3/4 ^" )u = = -/3(/3 - 1) + 3/4 - a 2 with allowed solution fi = a/I - a 2 + \ . At long distance 
(~lo% ~ x ~ fk) n = the- allowed solution has behavior u ~ x^ 1 ^" 2 cxp(— \x) which forces 7 = —1/2. 
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IV. THE TWO BODY DIRAC EQUATIONS OF CONSTRAINT DYNAMICS 

Dirac constructed a quantum wave equation from a first-order wave operator that is the matrix square-root of the 
corresponding Klein-Gordon operator [30] in order to treat a single relativistic spin-one-half particle, free or in an 
external field. The TBDE of constraint dynamics extend his construction to the system of two interacting relativistic 
spin-onc-half particles with quantum dynamics governed by a pair of compatible Dirac operators acting on a single 
16-component wave function. For an extensive review of this approach, see Refs. [4, 14, 17, 18] and works cited therein 
and [31]. We present below a brief review. 

Over thirty years ago, the relativistic constraint approach first successfully yielded a covariant yet canonical formu- 
lation of the relativistic two-body problem for two interacting spinlcss classical particles by applying a Hamiltonian 
approach introduced by Dirac [32] for handling systems with constraints. It accomplished this by introducing two 
constraints thereby reducing the number of degrees of freedom of the relativistic two-body problem to that of the 
corresponding nonrclativistic problem [33]- [36]. By this one covariantly eliminates the troublesome relative time and 
relative energy. The constraints used for this reduction are a pair of compatible generalized mass shell constraints for 
each of the two interacting spinless particles: 6 pf + mf + <f>i w 0. 

For the case of two relativistic spin-one-half particles interacting through four-vector and scalar potentials, the two 
compatible 16-component Dirac equations ([4, 14, 17, 18]) take the form 

Slip = 7 5 i (7i ' (Pi - M) + mi + Si)ip = 0, (26a) 
$2ip = 7 52 (7 2 ' (P2 - M) + m 2 + S 2 )iJj = 0, (26b) 

in terms of Si operators that in the free-particle limit become operator square roots of the Klein-Gordon operator. 

The relativistic four-vector potentials A% and scalar potentials Si are effective constituent potentials that in cither 
limit mi — > oo go over to the ordinary external vector and scalar potentials of the light-particle's one-body Dirac 
equation. The covariant spin-dependent terms in Af and Si are recoil terms whose general forms are nonperturbative 
consequences of the compatibility condition 

[5i,5 2 ]^ = 0. (27) 

This condition also requires that the potentials depend on the space-like interparticle separation only through the 
combination 

x^ = W + P»P v ){x 1 -X2) v (28) 

with no dependence on the relative time in the cm. frame. This separation variable is orthogonal to the total 
four-momentum 

P" (29) 

P is the time-like unit vector 

P^ = P^/w, (30) 
where w is the total cm. energy (the invariant rest mass), 

w 2 = -P\ 

so that in the cm. frame P = (1, 0) and x± = (0, r). The accompanying relative four-momentum canonically conjugate 
to x± is 

p^ = (e 2 P2 - £iP 2 )/w; where ei + e 2 = w, e± - e 2 = (ml - m 2 )/w. (31) 

The ej's are the invariant cm. energies of each of the (interacting) particles. Another consequence of the compatibility 
condition is that the relative momentum is constrained to be orthogonal to the total four-momentum 

P-pip = 0, (32) 



We use the metric r\ = (—1, 1, 1, 1). 
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thus providing the conjugate covariant control on the relative energy to that on the relative time provided by Eq. 
(28). One finds also that the vector and scalar potentials are defined in terms of two invariant functions S(r), A(r) 
in which r is the invariant 



r=^x\. (33) 

Those potentials have the general forms 

^ = i?(i4(r),px,P,«>,7i,72). Si = S i (S(r),A(r),p±,P,«',7i,72)- (34) 
The wave operators in Eqs. (26a) and (26b) operate on a single 16-component spinor 



V>= P (35) 




in which the ^ are four-component spinors. 

With compatibility ensured, this two-body formalism has many advantages over the traditional Bethe-Salpeter 
equation and its numerous three dimensional truncations. One is its simplicity. A Pauli reduction and scale trans- 
formation ([4, 14, 17, 18]) brings these equations to this covariant relativistic Schrodinger equation involving a four 
component spinor ip + , 



+ $ w {a u cT2,P±,A(r),S(r)) J V+ = & («#+, (36) 

resembling an ordinary Schrodinger equations with the interaction term & w including central-potential, Darwin, spin- 
orbit, spin-spin, and tensor terms. The interactions are completely local but depend explicitly on the invariant cm. 
total energy w — —P 2 . The usual invariant 

b 2 (w) = {w 4 - 2w 2 (m 2 1 + ml) + (to 2 - m 2 2 ) 2 )/Aw 2 , (37) 

plays the role of energy eigenvalue in this equation. This invariant is the cm. value of the square of the relative 
momentum expressed as a function of the invariant mass w. 

Note that in the limit in which one of the particles becomes very heavy, this Schrodinger equation turns into the 
one obtained by eliminating the lower component of the ordinary one-body Dirac equation in terms of the other 
component (when S(r) — 0, see Eq. (17) ). 



A. Relativistic Schrodinger Equation obtained from the TBDE 

In Appendix A we outline the steps needed to obtain the explicit cm. form of Eq. (36). That form is 

{p 2 + 2m w S + S 2 + 2e w A - A 2 + <S> D 

+ L • (cr 1 +<T 2 )(E> S o + <Ti-?<r 2 -?L • (trx+^^sor 

+ (Ti-a 2 ^ss + (3£rrf<T2-r-<x r er 2 )<I>T 

+ L • (eri-er 2 )$ SOZJ + iL-ai X cr 2 ^sox}^ + 

= b 2 ^ + . (38) 

The detailed form of the separate quasipotentials $i are also given in Appendix A. The subscripts of most of them 
are self explanatory 7 . 



7 The subscript on quasipotcntial <E>£> refers to Darwin. It consist of what are called Darwin terms, those that are the two-body analogue 
of the two terms to the left of the spin-orbit term in the one-body Pauli reduction given in Eq. (17), and ones related by canonical 
transformations to Darwin interactions [2, 37], momentum dependent terms arising from retardation effects. 
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After the eigenvalue b 2 of (38) is obtained, the invariant mass of the composite two-body system w can then be 
obtained by inverting Eq. ( 37). It is given explicitly by 

w = <Jb 2 +m 2 + ^Jb 2 +m 2 . (39) 

In the weak potential limit in which the potential is small compared with the masses and cm. energy, the effective 
potential <& terms of the relativistic Schrodingcr equation (38) contain singular potentials. Consider for example, the 
case of just a vector potential. The spin-spin term $55 includes through the part k(r) a piece |V 2 <7 (see Appendix 
A) that for weak potentials, i.e. for \A\ « w/2, is equal to V 2 A/(3w). (All other parts of $ss are negligible 
for weak potentials.) When placed as a correction to the potential in the nonrelativistic limit (so that it must be 
divided by 2fi) its contribution would be $55/2// = V 2 'A / (677117712) which, as discussed earlier in Eq. (9), would 
produce a singular delta function for a Coulomb- like potential. A similar analysis would show that for weak potentials 
$d/2^ = V 2 A/(8fi 2 ), a repulsive singular delta function potential for the Darwin term. For the spin-orbit term, 
$5o/2^ = A'/r(l/ (8/j 2 ) + 1/ (477117712)) which when combined with the spin-orbit eigenvalues produces either a 1/r 3 
attractive or repulsive singular potentials. One finds a similar behavior for the spin-orbit difference term. For the 
tensor terms one would find the combination ($sot — 3>t)/2/j = (A" — ^-)/(12m\m2) appearing in all angular 
momentum states producing attractive singular 1/r 3 potentials. (See radial equations ( B2)-(B4) below.) As we shall 
see below, when the strong potential structures are included, those type of singular potential structures disappear. 
However for the tensor terms we shall find below a remnant attractive singular potential (~ 7-~ 2 - 5 ) whose effect, 
however, is compensated by the tensor coupling. 

In the above equation (38) one has the strong potential form if one does not ignore the potential compared with the 
masses or the cm. energy in the expression for the $'s. In the following section we demonstrate that these strong 
potential forms and the unusual singularity structures in them lead to wave function solutions that are well-behaved 
and physically acceptable when the full couplings of the system are taken into account. We do this by examining the 
short distance behavior of these equations. 



V. THE BEHAVIOR OF SINGLET AND TRIPLET WAVE FUNCTIONS AT SMALL r IN QED 

In QED bound state problems, as has been shown in detail in [4] for A = —a/r the nonperturbative (numerical) 
treatment of the TBDE give the same spectral results as a perturbative treatment of the weak potential form (which 
in turn gave the same results as older standard methods). We note that this agreement was for a number of radial and 
orbitally excited states for equal and unequal mass as well as the ground state. In this section we examine analytically 
the structures of the four-component form of Eq. (38). 

For the equal mass singlet states, there is no spin-mixing and with A = —a/r, S = we obtain 



dr 2 



2e w a 



a 



vq — b v Q . 



(40) 



As was shown in [38] where this equation first appeared, it has an exact solution with the eigenvalue given by 



w = m 



\ 




l + a 2 



{L+\)*-c?-L-\ 



2m - ma 2 /An 2 - (ma 4 /2n 3 )[l/(2L + 1) - ll/32n] + 0(a e ), 



(41) 



or w — 2m — ma 2 /A — 21ma 4 /64 for the ground state. The small r wave radial function behaviors vq 



/■ 2 



v ' 1 ia ' 2 2 . Of course, its 



which implies the mildly singular but physically acceptable behavior of ip = Vo/r 
perturbative treatment gives the same results as the exact one through order a 4 . 

The ground state triplet cm. energy (excluding the annihilation contribution) perturbatively is [4] w — 2m — 
ma 2 1 A — ma 4 /192. This result, unlike the singlet case has not been obtained as an expansion of an exact analytic 
result. Unlike the one-body Dirac equation, the triplet states do not possess a known exact spectral solution for the 
TBDE. However, it has been verified that a nonperturbative (numerical) solution of the TBDE does produce a result 
that agrees with the perturbative evaluation [4] . 

We are more interested here in the behavior of the potentials in Eqs. (B3) and (B4) and the resultant wave functions 
at small r. Those, we show, can be determined analytically and provide a severe test for the strong potential forms 
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of the effective potentials. By using the effective potentials in those equations we find that at small r Eqs. (B3) and 
(B4) become (A = — a/r, S = 0) for a general J 



d?_ 

dr 2 



a (J+1)(J 

4(2.7+1) ) 2J(J + 1) 1 



+ 



2y/J(J + l) 
2J+ 1 



8r 



+ 



J + l 



a 
2rw 



2J+1 r 2 

- =0, 



a 
2rw 



u+ 



(42) 



and 



d?_ 

dr 2 



a 2 | J(J 



4(27+1)) | 2J(J+1) 



+ 



2y/J(J+l) 

2J+ 1 



" 1 J 


/-' 


- g r 2 ~ r 2 ^ 


/ 2rw 



u+ 



(2J+ l)r 2 V 2rw 
0. 



(43) 



Note that in Eq. (42), the potential for the u+ wave function contains an attractive term that is proportional to r -5 / 2 . 
Without the coupling to the U- wave function, Eq. (42) would lead to an attractive singular potential. However, 
the potential for the w_ wave function in Eq. (43) contains a repulsive term that is proportional to r~ 5//2 . The wave 
function u + is thus prevented from collapsing to the center due to the coupling of u + to m_ in Eq. ( 42). In fact, 
because of the coupling, the short distance behavior of u+ and u_ become proportional and given by 



u + (r) = r x 
u_(r) 



Jr x 



Ju+{r) 



y/J(J+l) y/J(J + l)' 



(44) 



which lead to an exact cancellation of the singular attractive r 5 / 2 terms in Eqs. (42) and (43) with their repulsive 
r~ 5 / 2 counterparts. The power index A with the correct physical behavior at the origin is 



A = v /J(J + l)- a 2 + i 



(45) 



so that at short distance the correct physical behavior of the wave functions is 



u+(r)=r (l/2+V^'+l)-« 2 ), 
Jr (l/2+VJ(i+l)-c 2 ) 

?i_(r) = 



The corresponding radial parts of the wave functions would be 

V> + (r)=7V-'(./+D-«M 7 

4>_(r) = 



(46) 



(47) 



Focusing on J = 1, both of these wave functions would show an unusual feature of a dip not only in the D state but 
also in the S state. This is in contrast to what occurs in the singlet case where there is the mildly singular behavior 
or in the nonrelativistic case where the behavior is flat. Also, whereas an upper bound of a = 1/2 is placed on the 
coupling in the singlet case, here the upper bound is \[2 for well-behaved wave functions. The unusual dips of the 5- 
and D-w&ve functions for QED also shows up in numerical solution of the QCD case, as will be seen in Figs. 3 and 4. 
For the special case of J = ( 3 Pq) the short distance behavior is 



r d 2 a 2 , 
{-^2-72-^ = 0, 



(48) 



just as it is for the 1 S'o case. 
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For the coupled 1 J,/, 3 Jj states, there will be spin mixing in the general case of unequal constituent masses. For 
equal mass or J = the equations decouple. However, in the limit of small r, the mixing term vanishes anyway, 
leaving us with the uncoupled short distance behavior for the wave equations (Bl) and (B2) for J > given by 

r d 2 J(J + 1) a\ 



d 2 J(J+1) a 2 1 

\- 

The radial wave functions behave as 



{-it + -^-^-^}vs=i = 0. (49) 



vs=o 



r (l/2 +A /(J+l/2) 2 -a 2 ) 



v s= i -> r (V2+V^(J+iP^). (50) 
The corresponding radial parts of the total wave functions would be 

^^..(-i/s+VU+Vs) 2 -^). 



fe ^ r H/2+\/W) ] (5i) 



and have acceptable behavior for J > when a is bounded by J +1/2 and y\7(J + l) respectively. Notice that unlike 
the tensor mixing case, there is no short distance connection in the spin mixing case between the wave function scales 
as in Eq. (44). For J = 0, only the first wave function is relevant and the wave function is well-behaved at short 
distances for a bounded by 1/2. 

We can summarize the short-distance behavior of the wave function for different spin and angular momentum states. 
For the coupled or uncoupled ( 1 Jj, 3 Jj) states, the corresponding potentials arc not singular. In the tensor coupling 
case, there are both attractive and repulsive singular potentials. However, the effect of the coupled equations is to 
produce eventually well-behaved wave functions. Thus all wave functions at short-distance are well-behaved for all 
spin and angular momentum states for appropriate a bounds. 

Next we examine the singularity structure of the TBDE in the case of QCD bound states. The short distance 
behaviors (which we attribute to the invariant A(r)) of QED and QCD are well known to have a crucial distinction. 
Renormalization group arguments show that in QED the asymptotic behavior displays a singularity structure 8 at 
finite (though very large) energy that invalidates perturbation theory [40] . In contrast similar renormalization group 
arguments show for QCD the asymptotic behavior displays a structure 9 that at high energy strongly validates per- 
turbation theory (asymptotic freedom) [40] . For this reason we have not included a running coupling constant above 
in the QED application. 



VI. SIMPLE QCD BASED MODEL FOR QUARK-ANTI-QUARK BOUND STATES. 

In [17] a fully relativistic calculation of the meson spectrum was made using the above TBDE with the invariants 
A(r) and S(r) determined from a relativistic extension of the nonrelativistic Adler-Piran [19] static quark potential. 
That paper investigated how well the relativistic constraint approach performs in comparison with selected alternatives 
including those of Godfrey and Isgur [41] when used to produce a single fit of experimental results over the whole 
meson spectrum. The authors of [17] found that the fit provided by the two-body Dirac model for the entire meson 
spectrum competes with the best fits to partial spectra provided by the others and does so with a smaller number of 
interaction functions. Furthermore this is done without additional cutoff parameters necessary to make some of the 
other approaches numerically tractable. 

In this section we examine the spectral results produced by a relativistic version of a simpler model motivated 
by the static quark potential of Richardson [23]. It is not our purpose here to improve upon the results of [17] but 
rather to take advantage of the simpler structure of the model to explore certain details of the effective potential and 
wave function behaviors generated by the TBDE along the lines discussed in the above sections for QED. This will 
also allow us to more readily make changes in the invariants A(r) and S(r) determined from this simpler model to 
facilitate investigations on quarkonium stability in the presence of a quark-gluon plasma. 



8 The running coupling constant in spinor QED is given by [39] an(q 2 ) — ► o/j(m 2 )/(l — a^(m 2 ) ln(— q 2 /m 2 )/3n). 

9 The running coupling constant in QCD is given by [39] l/a s (<j 2 ) — > l/a s ()i 2 ) + (33 - 2n f)\n{-q 2 / fi 2 ) /Y2ir . 
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A. The Model and Spectral Results 



For our potential model we consider here a slightly more generalized form of Eq. (8), 



V(r) 



8vrA 2 



16tt 



27 27r\n(Ke 2 + B/(Ar) 2 )' 



(52) 



This has the feature of giving a variable parameter K to the long distance Coulomb behavior as well as (through B) 
an effective QCD parameter A distinct from the one placed in the linear potential. In our QCD spectral work we 
assign the confining piece of this potential to the invariant S that controls the scalar potential and the Coulomb-like 
piece to the invariant A that controls the vector potential (see Eqs. (A2 ),(A5),(A6), and (All)). Thus with r = \l x\ 
the equations 



S(r) = 



8ttA 2 



27 



A(r) = 



16?r 



eie 2 



27rln(ATe 2 + 5/(Ar) 2 ) Airr 



(53) 



together with the relativistic Schrodinger equation (38) of our TBDE (A12) define the covariant formalism for our 
QCD spectral work (ei,e2 are the respective electric charges of the quark and anti-quark). With the parameter 
values listed in Table I we obtain the following spectral results shown in Table II and Figure 1. The mass of the u 
and d quarks are only 55 MeV. We observe that the masses of the low-lying mesons are well reproduced in the TBDE 
treatment with a minimum number of parameters. The singlet-triplet splittings of ir-p and rj c -J /ip are well reproduced. 
The agreement of the experimental masses with theory is not as impressive as that using the Adler-Piran potential. 
(Note that as with the Adler-Piran potential, the ground state singlet/triplet splitting for the charmonium system 
appears too large, while that between the ground and first excited states appears too small. The same problem also 
occurs for the bottomonium system with the recently observed r\ h [42].) The simplified potential has the advantage 
of simplicity and ease of adoptive modification that can be useful for applications of the TBDE to other quarkonium 
problems. 
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FIG. 1: (Color online) Comparison of experimental and theoretical masses obtained with the Two-Body Dirac Equations. 



We next consider some interesting and unexpected behaviors of the solutions of the equations in the next few 
subsections. 
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A = 


= 0.4218 Gev 


B = 


= 0.05081 


K -- 


= 4.198 


m u 


= 0.0557 GeV 


m d 


= 0.0553 GeV 


m 3 


= 0.2499 GeV 


m c 


= 1.476 GeV 


m b 


= 4.844 GeV 



TABLE I: Quark Model Parameter Values 



|\/Tpcnn 

IVlCoUll 


Exp. (GeV) Theory (GeV) Exp.-Theory(GeV) 


7T . Vb(l 1 D() 


0.14U 


U.159 


-o.oiy 




0.775 


0.792 


-0.017 


K : sii 1 bo 


U.494 


n a no 


U.001 


A : sa 1 on 


0.498 


0.488 


0.010 


K : su 1 hi 


0.892 


0.903 


-0.011 


K* : sd 1 3 5i 


0.896 


0.901 


-0.005 


: s s 1 3 5i 


1.019 


1.025 


-0.006 


r-.fl -tin 


1.865 


1.840 


0.025 


D + : erf l^o 


1.870 


1.845 


0.025 


V : cu 1 6i 


2.010 


1.981 


0.029 


D* + : cc/ 1 3 Si 


2.007 


1.979 


0.028 


_D S : cs 1 So 


1.968 


1.965 


0.003 


L) s . cs 1 oi 


119 


110 


u.uuu 


- 1 1 c 

Tj c \ CC L JO 


2.980 


2.978 


0.002 


J/i>(lS) : ccl a Si 


3.097 


3.140 


-0.043 


tp(2S) : cc 2 3 Si 


3.686 


3.689 


-0.003 


hi : cc 1 1 Pi 


3.526 


3.522 


0.004 


X : cc 1 3 P 


3.415 


3.436 


-0.021 


Xi : cc 1 3 Pi 


3.511 


3.515 


-0.004 


X 2 = cc 1 3 P 2 


3.556 


3.541 


0.015 


r; c : cc 2 So 


3.638 


3.591 


0.047 


V)(1D) : cc l 3 Di 


3.773 


3.804 


-0.031 


B - : bul'So 


5.279 


5.249 


0.030 


B° : bd 1 1 S 


5.280 


5.248 


0.032 


B*° : bu l 3 Si 


5.325 


5.299 


0.026 


B s : bsl 1 S () 


5.366 


5.360 


0.006 


B* 3 : bs 1 S S! 


5.413 


5.420 


-0.007 


B~ be 1 1 So 


6.276 


6.276 


0.000 


V b ■ bbl'So 


9.389 


9.345 


0.044 


T(15) : 6b l 3 Si 


9.460 


9.484 


-0.024 


T(2S) : bb 2 3 Si 


10.023 


10.033 


-0.010 


T(35) : bb 3 3 Si 


10.355 


10.360 


-0.005 



TABLE II: Selected Portions of Meson Spectrum 



B. Behaviors of singlet and triplet solution to bound state equations 

1. Detailed Analysis of the Pion and Rho Bound States. 

One of the most unusual features of the spectral results of the TBDE is the fairly accurate production of the large 
-K-p mass splitting. First we point out that unlike most other potential models, the up and down quark masses do 
not take on the values typically seen of about 300 MeV in most other potential models (see e.g. [22]). Just as with 
the more detailed Adler-Piran potential [19], the quark masses that give the best fit in our model here are on the 
order of 50 — 60 MeV. It is of interest to see how the small pion mass comes about. In earlier work [17], it was shown 
numerically that as the quark masses tend to zero, so also does the bound state mass. Thus our analysis here bears 
on the dynamics of chiral symmetry breaking. The bound state equation for S states has the general form of Eq. (Bl) 
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and in this case reduces to 

d 2 

{--^2 + 2m ^S + S 2 + 2e w A -A 2 + <P D -3<S>ss}v 

d 2 d 2 

= {-71 + *l v o = {-77 + 2m ^ s + s2 + 2e ^ A -A 2 + <S> SD } 
ar dr z 

d 2 

= {—-^ + <S>si + ®SD} = b 2 v . (54) 

Table III below lists the expectation values of the various parts (units in GeV 2 ) . Even though the quark masses make 
up a substantial portion of the pion mass it would be a gross mis-statement to say that the pion bound state is nearly 
nonrelativistic. The constituent kinetic portion (— -4-%) and potential portion ($) are huge compared with the rest 
mass squared (0.00302 GeV 2 ) but nearly cancel, leaving a small b 2 which corresponds to a pion mass of 0.159 GeV. 
The pion as a quark- antiquark system has large kinetic energies and potential energies that counterbalance each other. 
When one looks at the potential energy ($) = —0.8475 GeV 2 from various contributions, one notes that the spin- 
independent contribution ($>si) = (2m w S) + (S 2 ) + (2e w A) - (A 2 ) = -0.3832 GeV 2 while the combined Darwin and 
spin-dependent contributions are ($sd) = (&d — 3$ss) =-0.4643 GeV 2 . They are of the same order of magnitude. 
Among the two huge contributions to (&sd), the Darwin term (<&d) and (— 3$ss) nearly cancel one another. The 
small mass of the pion arises in no small measure from delicate cancellations of these very large contributions from 
the Darwin and spin-spin interactions, as well as the balance of kinetic energy and potential energies. Such does not 
occur with the other mesons. (See for example the rho meson below). As seen in Fig. 2 there is nothing unusual in 
the pion wave function in constraint dynamics. 




FIG. 2: (Color online) Wave functions of n and r\ c mesons 



By contrast now we present a similar table for the p meson. The coupled equations are 



d 2 

+ 2m w S + S 2 + 2e w A - A 2 + $ D + $ss}u+ 

/o 

+ — {3$ T - 6$ SO t}u- 
j- d 2 

= i"^2 + + 

= b 2 u +1 (55) 



and 



15 



($> 


0.8508 
-0.8475 


(2m w S) 


0.0103 


<^> 


0.0942 




-0.0598 


<-^> 


-0.4279 


<$s/> 


-0.3832 


<<M 


-3.804 


<-3$ss) 


3.340 




-0.4643 


<^> 


0.0033 



TABLE III: Expectation values of various terms in (54) (in GeV 2 ) that contribute to the pion eigenvalue equation (54). 



d 2 6 

{--To + ^ + 2m wS + S 2 + 2e w A - A 2 + <P D - 6$ so + &SS - 2$ T + 2$ SO t}w- 
dr z r z 

+ — {3$ T }u+ 

r d2 6 * 1 
= {-^2 +^2+ $ — }«- + *-+«+ 

= fe 2 u_. (56) 
In terms of expectation values we have 



r°° f (f 

b 2 = J rf r[w + (r) -— + $++) u + (r) +u + (r)$ + _-fi_(r) 



( d 2 6 \ 
+ u -( r ) +^2+ $ — J«-W+«-W*-+«+W] (57) 

Table IV below lists the expectation values of the various parts (units in GeV 2 ). The quark masses make up only 
a small portion of the rho mass. As with the pion, the rho is highly relativistic, but unlike the pion, its relativistic 

nature is not hidden in large cancellations. The constituent kinetic portion (—^3)++, {—-§pz + \) , and potential 

portions ($)++, (<&)_| ($) |_, ($)— are huge compared with the rest mass squared (0.00302 GeV 2 ). There is some 

cancellation but not nearly to the extent that occurs in the pion. There is left a significant b 2 =0.1411 GeV 2 which 
corresponds to a computed rho mass of 0.796 GeV. Note that by itself, the 5-wave portion (kinetic plus potential) is 
negative and large compared to the rest mass. Thus the positive D-wave portion is crucial to bring the rho mass in line 
with the observed value. The two spin-independent contributions are ($>si)u — (2m w S)u + (S 2 )u + (2e w A)u — {A 2 )a 

where ii = ++ and . The magnitudes of ($<?/)++ = -0.1680 GeV 2 and ($<?/)— = 0.01461 GeV 2 differ by an 

order of magnitude, as do the diagonal Darwin and spin-dependent contributions ($sd)++ = ($d + ®ss)++ =-0.3154 
GeV 2 and ($ S d)— = (<&d + *S5 - 2$t - 6$ so + 2$ SO t)-- = 0.03461 GeV 2 . The two sets differ by roughly a factor 
of 2 and arc of the same sign. The off diagonal tensor terms ($)+_ = 2 -j^(3$t - 6$sot)+-, ($)-+ = 3^(3$t)-+ 
are both quite large but of opposite sign so their overall effects almost cancel. 



2. Detailed analysis of the vector meson potentials and wave functions. 



As may be anticipated from the dip behavior displayed in Eq. (47) for the case of QED interactions we anticipate 
that a similar structure may appear for the QCD wave functions. This is borne out by Figs. 3 and 4 for the S and 
D state contributions to the J/ip and the p mesons. In the equations below the analytic origin of the dip behavior 
is shown from the behaviors of the quasipotential contributions. 

In contrast to the short distance behavior displayed in Eqs. (42) and (43) in the QED case, the short distance 
behavior for our coupled QCD equations, corresponding to mesons such as the J/ip and p mesons is 
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(-;£)++,(--£,. + 4) — 


0.3085 0.2812 


($)++, (*)_. 


-0.4835 0.04923 


<$)+_, {$)_+ 


0.3090 -0.3088 


(2m w S)++, (2m w S)— 


0.00263 0.000571 


X" /++> W /— - 


1631 04457 


/Op 4\ , , /9p A) 


-0 2091 -0 02109 


\ / H — r 3 \ -* 1 / 


-0.1247 -.00944 


($SI>++, ($SI> — 


-0.1680 0.01461 


($ D )++, ($d) — 


-0.2790 -0.04360 


($Ss}++, ($Ss) — 


-0.03637 -0.00947 


(-2*T>- 


-0.04172 


-6(*so)- 


0.1133 


2(<E> SO t)-- 


0.02153 


\*SI)/++i \y?SD) — 




^<3$ T > + _,2^(3$ T )- + 


-0.3088 -0.3088 


^2(-6<I> SOT > + _ 


0.6177 


<$>+_, <<&)_+ 


0.3090 -0.3088 


<^> ++ , <&*)__ 


-0.1750 0.3158 




0.3090 -0.3088 


(6 a )T taI=0.1411 





TABLE IV: Expectation values of various terms in Eqs. (55) and (56) (in GeV 2 ) that contribute to the rho eigenvalue equations 
(55 and (56)). 
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r(fm) 

FIG. 3: (Color online) Wave functions for the J/ip meson. The insert in the upper right corner gives an expanded view of the 
wave functions near the origin. 
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27u>rln(Ar) 6r 2 



2y/2 
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27w ln(Ar-) 8r 2 
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(58) 
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FIG. 4: Wave functions for the p meson. The insert in the upper right corner gives an expanded view of the wave functions 
near the origin. 



and 



d 2 6_ 

dr 2 r 2 3r 2 

= 0. 

Let us assume a short distance behavior of 



4tt 49 


U- 


2y/2 




27wr ln(Ar) 12r 2 


f 3 


r2 V 



4tt 



r 2 V 27wrln(Ar) 8r 2 



/(r) _ u+(r) 



U-(r) = 



V2 



(59) 



(60) 



Then just as in the case of QED, the more singular terms (here ~ (I / r 2 ) —Air /27wr ln(Ar)) cancel among themselves 
and we are left with 



and 



which have the identical behavior of 



" d 2 


11 " 


m- 


2 


1 


f(r) 




dr 2 


6r 2 


h 3 


"Ir 2 


= 0, 


' d 2 


23 


f(r) 


4 


1 


f(r) 




dr 2 


12r 2 


+ 3 




= 



dr 2 4r 2 



f(r) = 0, 



fM ~ r (V2+V2) 

The corresponding radial parts of the wave functions would be 

V+(r)=r^-5, 



(61) 

(62) 

(63) 
(64) 



V2 ' 



(65) 
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compared to the corresponding radial parts of the QED wave functions of 



(66) 



V2 



Numerically, one finds results that approach those of Eq. (65) in terms of the S/D ratio and power behavior. 



We have shown how the TBDE of constraint dynamics handles the problem of effective potentials that are singular 
in the weak potential limit. The most noteworthy feature is that different portions of the quasipotential $ contribute 
to the perturbative and nonperturbative treatment of spectral effects. The results we found for the coupled 3 S\- 3 Di 
system demonstrates that this effect extends not only to different terms of a given equation, but also bridges the 
divide between the various coupled components of the wave function and effective potential. We find unexpectedly, 
that the behavior of the J = 1 , S and D waves (and more generally for an arbitrary J) are simply proportional very 
near the origin, with a common power-law behavior. 

We have introduced a new QCD potential that has many of the features of the Adler-Piran and Richardson 
potentials but with a much simpler parametrization. Although not giving results (rms deviation=21 MeV) as good 
as the former [17] (rms deviation=14 MeV), it nevertheless yields a meson mass spectrum that agrees reasonably well 
with experiment. We examine in particular how various contributions to the pion eigenvalue equation can lead to a 
pion of a small rest mass. The detailed treatment of the pion shows a unique feature of it relativistic behavior, namely 
a behavior that superficially appears nonrelativistic but on further analysis displays extremely large relativistic and 
nearly canceling contributions to and from the potential. This gives some insight into how the potential model leads 
to a small pion mass for small quark mass, an important consequence of spontaneous symmetry breaking. By contrast 
our Table IV shows why mesons with other quantum numbers, would not display such a small quark mass behavior 
for the bound state energy. 

Another unexpected behavior for the QCD wave functions and effective potentials is that at short distance, the 
J = 1, S and D waves for the J /ip mesons are not significantly different from that which appears in QED, in spite of 
the asymptotic freedom behavior that occurs in QCD. It may be worthwhile to investigate what observable features 
of these bound states (both in QED and QCD) may reflect these unexpected connections between the two coupled 
wave functions. 

In order that our primary bound relativistic Schrodingcr equation (38 ) and the radial equations (B1)-(B4) have 
an additional ease of use for others who may wish to apply them, we have added Appendix C to detail the relation 
between the $'s that appear in (A18) and the invariants A(r) and S(r) and their derivatives. 

In addition to applications to meson spectroscopy, our work has implications for studies related to quarkonia at 
finite temperatures and in a quark-gluon plasma [43-46]. Within the potential model, the quark drip line at which 
quarkonia begins to be unbound has been estimated using non- relativistic quark models [47, 48]. A more definitive 
investigation of the stability of quarkonium states, especially those with light quarks, necessitates the use of the 
relativistic formalism developed here. The analytic structure of the A and S invariants given in Eq. ( 53) lend 
themselves to practical modeling at finite temperatures, while giving an adequate account of the zero temperature 
limit. The relativistic one developed here provides a more realistic model than the most often used Cornell-type 
linear-plus color-Coulomb for the case of T = 0. It will be of interest to investigate the composite properties of the 
plasma and to determine the region of temperatures in which different quarkonia become unbound using a temperature 
dependent extension of the simple quark model presented in this paper. An important problem for future work in 
relation to the present paper is to determine how best to apportion these temperature dependent potentials between 
the A(r) and S(r) invariants. 
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VII. CONCLUSION AND FURTHER REMARKS 



APPENDIX A: RELATIVISTIC SCHRODINGER EQUATION DETAILS 



Here we present an outline of some details of Eq. (38) given in full elsewhere (see [4, 14, 17, 18, 49, 50] and works 
cited therein and [31]). For classical [2] or quantum field theories [3] for separate scalar and vector interactions one 
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can show that the spin independent part of the quasipotential involves the difference of squares of the invariant 
mass and energy potentials (Mj and E t respectively) 

Mf = m 2 + 2m w S + S 2 ; E 2 = e 2 - 2e w A + A 2 , (Al) 
M 2 - E 2 = 2m w S + S 2 + 2e w A - A 2 - b 2 (w). (A2) 

where 

_miTO2 _ (w 2 -m\- ml) 

w 2w 

are respectively the relativistic reduced mass and energy of the fictitious particle of relative motion introduced by 
Todorov [1, 12] and satisfy the effective one-body Einstein condition 

e 2 w -m 2 w = b 2 (w). (A4) 

Eqs. (26a) and (26b) contain an important hidden hyperbolic structure [51]. To reveal and employ it one introduces 
two independent invariant functions L(x±) and Q(x±), in terms of which the invariant mass and energy potentials 
take the forms: 

Mi — mi coshL(S,A) + m 2 smhL(S, A), 

M 2 = m 2 coshL(S,A) + m x sinhL(S,A), (A5) 
E 1 =e 1 coahG(A) - e 2 sinh^(A), 

E 2 = e 2 coshG(A) - si sinh£(A). (A6) 

Strictly speaking, the forms in Eq. (Al) are for scalar and time-like vector interactions. Eq. (38) below involves 
combined scalar and electromagnetic-like vector interactions (this amounts to working in the Feynman gauge with 
the simplest relation between space- and time-like parts, see [8, 17])). In that case the mass and energy potentials in 
place of Eq. (Al) are respectively 

M 2 = mj + atv(2G)(2m w S+S 2 ), (A7) 

so that 



exp(L(5 A)) = V^i + cx P( 2 g) (2^5+S 2 ) + y/m 2 + cxp(2g)(2m^+S 2 ) 

mi + m 2 

and 

£ 4 2 = cxp(2g)( £t -A) 2 , (A9) 

with 

exp(W)) =JT*am 55 ° 2 - (A10) 

In terms of Q and the constituent momenta pi and p 2l the individual four-vector potentials take the forms[17] 

1 

Ai = [1 -cosh(^)]pi + sinh(£)p 2 - -(dexpG ■ j 2 )-y 2 , 

A 2 = [1 - cosh(G)}p 2 + sinh(g)pi + l -{dexpQ -7ihi- (AH) 
In terms of the three sets of invariants (A7)-(A10) the coupled TBDE (26a) and (26b) then take the form 
Snl> = ( - GftEi • V 2 + Ei(3i l5 i + Mi 761 - G l -Y, 2 ■ 3{GPi + £/3 2 )75i7 52 )^ = 0, 
S 2 V = (G/3 2 E 2 • Vi + E 2 p 2l52 + M 2752 + G l -^i ■ 8(GP 2 + L/3 1 )75i7 52 )^ - 0, (A12) 



20 



in which 



Vi=p - ^i-dQ^i, (A13) 



depending on gamma matrices with standard block forms (see Eq. (2.28) in [17] for their explicit forms) and where 

= 7 5i &7±i- (A14) 

The Klein-Gordon like potential energy terms appearing in the Pauli form ( 38) of Eq. (A12) are as in Eq. (A2). 

To obtain the simple Pauli form of Eq. (36) involves steps analogous to those used in Eqs. (11)-(13) but with the 
combinations <j>± = ip 1 ± ip i and \± = "02 ^ ^3 instead of the the individual ip i . This allows the Pauli forms to reduce 
to 4 uncoupled 4 component relativistic Schrodinger equations [17, 18, 50, 52, 53]. Working in the cm. frame in 
which P = (1, 0) and f = (0, f) and then further defining four component wave functions ip±,r]± related to the above 
by [53] 

<j>± — cxp(F + K(Ti-r(T2-r)tp ± = (exp F)(cosh K + sinh K ovf <x 2 f )ip±, 

X± — cxp(F + K(Ti-r(T2-r)v± = ( ex P F) (cosh K + sinh.fi: <x r fer 2 -f)?7 ± . (A15) 

in which 

1 V 

F = -i og g, 

2 e 2 mi+£im 2 
V=E 2 Mi + E 1 M 2 , 

k = ^ (A16) 

will yield equations which in the limit when one mass becomes extremely large reduce to the Schrodinger or Pauli- 
forms discussed in Section III. In analogy to what occurred there the decoupled form of the Schr odinger equation for 
ip, has the convenient property that the coefficients of the first order relative momentum terms vanish. 

Using the results in [53] we obtain for the general case of unequal masses the relativistic Schrodinger equation (38) 
that is a detailed cm. form of Eq. (36). In that equations we have introduced the abbreviations 

2F' (cosh2.fi: - 1) , 2 ~ 2K' sinh2.fi: „ 2 ^ 
$ D = 1 L + F + K' 2 H V 2 F 

r r 

2(cosh2if- 1) 

F' F'(cosh2K- 1) (cosh2.fi: -1) if' sinh2.fi: 
$ so — 1 5 

®sod = {I' cosh 2K - q sinh 2K) , 
®sox = (q' cosh 2K- /'sinh 2K), 

2K'sinh2K 2F '(cosh 2K - 1) 2(cosh2if-l) 

2F'K' V 2 K fK . 

+ — 3"' (A17) 

1, , s iF'sm\\2K F'(cosh2K - 1) , , lTsinh2.fi: 3if'(cosh 2K - 1) 

*t = -=[n(r) + + — '- + 2F'K' ^ '- 

3 r r r r 

2 3sinh2.fi: (cosh2.fi: -1) 

~~ V K H ^2 1 ^2 J' 

.cosh2.fi:- 1 sinh2.fi: K' F' smh2K , k N 

$SOT = ~K' + + . A18 

Hr^^iK^-'-^jn^ 
»(r) = i[V 2 if - iv 2 (? + 3(G ' 2r 2K,) + F'(G' 2K% 

m(r) = -±V 2 g+'±g' 2 + g'F'-K' 2 , (A19) 



in which 



21 



and 

l_E 2 M 2 -E lMl 
V ' 2rE 2 M 1 + E 1 M 2 ' 
1 E 1 M 2 -E 2 M l fT 
q{r) = Yr E 2 M 1+El M 2 {L - G) - (A20) 

(The prime symbol stands for d/dr). For L = J states, the hyperbolic terms cancel and the spin-orbit difference 
terms in general produce spin mixing except for equal masses or J = 0. For utility of use we have listed in Appendix 
C the explicit forms that appear in the above <&'s in terms of the general invariant potentials A(r) and S(r). The 
radial components of Eq. (38) are given in Appendix B. 

APPENDIX B: RADIAL EQUATIONS 

The following are radial eigenvalue equations corresponding to Eq. (38 ) [53]. For a general singlet 1 Jj wave function 
v coupled to a general triplet 3 Jj wave function v\, the wave equation 

{"^2 + ~ + 2m ™ S + s2 + 2e ^A -A 2 + $d-3$ss}«o 

+ 2^J{J + 1)($sod - <S>sox)vi 

= b 2 v , (Bl) 



is coupled to 

{-TT + J(J 2 +1) + 2m w 5 + S 2 + 2e w A - A 2 + <& D 
dr A r A 

- 2$ so + *ss + 2$ T - 2$sot}vi + 2y/j(J + T)($sod + $sox)v 

= b 2 v x . (B2) 
For a general S = 1, J = L + l wave function u + coupled to general S = 1, J = L — 1 wave function m_ the equation 

{ — 7-7 + J< ^ J 7 ~ + 2m™ S + S 2 + 2e w A - A 2 + $ D 

+ 2(J - l)* so + + ^ j + f (^sor - *t)}m+ 

2V%/TT) r 
+ 2 J + 1 {3$r " 2(J + 2)$sot}u - 
(,2 



is coupled to 



6 2 it+, (B3) 



d 2 , (J+l)(J + 2) , 



{-3^ + ^ + 2m w S + S 2 + 2e w A - A 2 + $ 

dr z r A 

2(J + 2), 

- 2(J + 2)$so + *S5 + 2J+l (®sot - $t)}u- 

+ 2 ^2J J + + l l) { ^ T + 2(J " 1 )*sor}«+ 
(,2 



= ^u_. (B4) 



22 



APPENDIX C: EXPLICIT EXPRESSIONS FOR TERMS IN THE RELATIVISTIC SCHRODINGER 

EQUATION (38) FROM A(r) AND S(r) 

Given the functions A(r) and S(r) for the interaction, users of the relativistic Schrodinger equation (38) will find it 
convenient to have an explicit expression in an order that would be useful for programing the terms in the associated 
equation (A18). We use the definitions given in Eqs. (A16), (A19), and (A7 )-( A10). In order that the terms in Eq. 
(A18) be reduced to expressions involving just A(r), S(r) and their derivatives, we list the following formulae: 

pi = {L' - Q'){E 2 M 2 + EiMj) , 
2(E 2 M 1 + E X M 2 ) 



Mi 
L' 




2m w S + S 2 
{w-2A)/w ' 

(2m w S + S 2 )A' 



(w - 2Af 

■ = —T7- (CI) 



Also needed arc 



cosh2if 



1 / (ei + £2) (Mi + M 2 ) (mi + m 2 )(E 1 + E 2 ) 



2 \{mi +m 2 )(Ei + E 2 ) (ei + e 2 )(M 1 + M 2 ) 



■ h o K = 1 / (£i+£ 2 )(M 1+ M 2 ) (m 1+m2 )(E 1+ E 2 y 
bm 2\{m 1 +m 2 ){E 1 +E 2 ) ( £l + e 2 ){M x + M 2 ) > ' 



and 



Q' + L' 

2 ' 

(V 2 £~ V 2 g)( J E 2 M 2 + g 1 Mi) _ _ 2 {m\ -ml) 2 _ 2( 

2(E 2 M 1 +E 1 M 2 ) ( y) 2{E 2 M l +E l M 2 f 

-L' 2 (M 2 + M|) 



MiM 2 



| w / V' i 5'(m tt , + 5) + 5' ,2 | 4SK + ^' + (2m w S + S 2 )V 2 A | 4(2771^5 + S 12 )^ 2 



V 2 L 

^MiM 2 V w-2A ' (w-2A) 2 ' (w - 2Af 

V * S = -^lA +2e °- < C3 > 

The expressions for k(r),m(r), and n(r) that appear in Eq. (A18) are given in Eq. (A19). They can be evaluated 
using the above expressions plus 

v3 ^VWL (C4) 

The only remaining parts of Eq. (A18) that need expressing are those for I' and q' . Using Eq. (A16) they can be 
obtained in terms of the above formulae. 
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